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Abstract 

In this article, we show that a stochastic form of scale-dependent halo bias arises in multi- 

source inflationary models, where multiple fields determine the initial curvature perturbation. 
^T) We derive this effect for general non-Gaussian initial conditions and study various examples, 

such as curvaton models and quasi-single field inflation. We present a general formula for both 
p\J the stochastic and the non-stochastic parts of the halo bias, in terms of the A'^-point cumulants 

0^ of the curvature perturbation at the end of inflation. At lowest order, the stochasticity arises if 

the collapsed limit of the four-point function is boosted relative to the square of the three-point 
I function in the squeezed limit. We derive all our results in two ways, using the barrier crossing 

J> formalism and the peak-background split method. In a companion paper [1], we prove that these 

^ two approaches are mathematically equivalent. 
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1 Introduction 



A central goal of modern cosmology is to uncover the physics that generated the primordial 
density perturbations and thereby seeded the large-scale structures (LSS) we see around us. 
The coherent nature of the cosmic microwave background (CMB) anisotropics suggests that the 
fluctuations were created at very early times, possibly during a period of inflation [2]. 

One of the few observational probes that allows us access to the physics of that epoch is pri- 
mordial non-Gaussianity [3]. At present, the best constraints on non-Gaussianity are coming from 
the CMB (e.g. [4]), but LSS is emerging as a promising complementary observable (e.g. [5, 6]). 
Historically, the usefulness of LSS as a tool for early universe cosmology has been viewed with 
some suspicion, since non-linear evolution can itself produce significant non-Gaussianity even if 
the initial conditions were perfectly Gaussian. Disentangling any primordial non-Gaussianity 
from these late time effects always seemed like a messy business. This attitude has changed 
somewhat when it was discovered that non-Gaussian initial conditions lead to a scale-dependent 
clustering of galaxies on large scales [7, 8]. In particular, it was shown that non-linear mode 
coupling induces a modulation of the local short-scale power as{x) by the long- wavelength grav- 
itational potential $(a;). This results in a biasing of halos (or galaxies) that is proportional to ^ 
rather than the dark matter density 6 (or V^$). Crucially, the appearance of $ rather than 6 
in the halo bias implies a specific form of scale-dependence that cannot be created dynamically 
(i.e. by late time processes). This is the main reason that halo bias is such a robust probe of the 
initial conditions. 




Figure 1: The squeezed limit of the three-point function, ki — > 0, gives the dominant contribution to the 
scale-dependent halo bias. A stochastic form of scale-dependent halo bias arises if the four-point function is 
large in the collapsed limit, ki2 = |fci + ^2! ^ 0. 



In this paper, we study stochastic halo bias on large scales. The term 'stochastic' here refers 
to the fact that the halo over-density is not 100% correlated to the matter over-density on large 
scales, i.e. the halo-halo power spectrum Phh(^) is boosted relative to the matter-halo power 
spectrum P^hik)- Formally, this means that 

Pthik) > h'{k)P^^{k) + — , (LI) 

nh 

where h{k) = -Pmh(^)/-fmm(^) is the halo bias, and is the halo number density. Large-scale 
stochastic bias arises in non-Gaussian models when the small-scale power a^{x) varies from point 
to point, but in a way that isn't completely correlated with the local value of $(a;) and its deriva- 
tives. This is most easily demonstrated in models with multiple fields, where the small-scale power 
may depend on fields that do not contribute to the (linearized) gravitational potential. Our goal 
in this paper is to provide an understanding of the origin of stochastic bias in a model-independent 
way. In the absence of significant isocurvature perturbations, all the relevant information must 
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be encoded in the correlation functions of gravitational potential It will be useful to define 

f _ 1 eg^(fci,fc2,fc3) 

^-=4^^o P[P, ' ^'-^^ 

where • • • '3>fejv)c = (27r)^^^^^(fci, • ••,/!;„) (5d(A;iH h^n) and Pj = C^Hh)- This parametrizes 

the amplitude of the three-point function in the squeezed limit, ki <^ min{A;2, Ajs}, and the am- 
plitude of the four-point function in the collapsed limit, ki2 = |fci + A;2| ^ minj/cj}. As we will 
show, stochastic bias arises if the 'collapsed four-point function' is not equal to the square of 
the 'squeezed three-point function', i.e. if Tnl 7^ (f /nl)^- There exists a well-known theoretical 
constraint on the relative size of Tnl and (|/nl)^- If only a single field (which may or may not 
be the inflaton) generates the primordial curvature perturbation and its non-Gaussianity, then 
Tnl = (|/nl)^ [9] and the biasing is non-stochastic. On the other hand, if multiple coupled fields 
generate the non-Gaussianity, then be larger^ than (|/nl)^ [16, 17] and the biasing will 

be stochastic. We will discuss classes of inflationary theories that predict precisely this kind of 
observational signature [18-20]. This provides the opportunity of using scale-dependent stochas- 
tic bias^ as a probe of any early universe physics associated with a boosted collapsed four-point 
function — ^just like the non-stochastic scale-dependent bias is a powerful probe of the squeezed 
three-point function. 

More generally, we find that the large-scale non-stochastic bias can be written as a sum 
over A^-point functions ^^\ki, • • •, k^) evaluated in the squeezed limit ki <C min{/c2, • • •, k^}-^ 
The stochastic bias, on the other hand, involves a double sum over (M + A^)-point functions 
C^^~^'^\ki, • • - jkM+w) evaluated in the collapsed limit jfci + • • • + fcA/| <^ minjfcj}. Stochastic 
bias arises if any collapsed (M + A^)-point function is boosted relative to the product of the 
corresponding squeezed (M + l)-point and (A^ + l)-point functions. In all physically interesting 
cases that we are aware of, this effect is due to the collapsed four-point function being boosted 
relative to the square of the three-point function (i.e. the case M = A^ = 2). Therefore, we 
will generally interpret stochastic bias as a probe of the collapsed four-point function. The main 
result of this paper is a general pair of formulas, eqs. (3.26) and (3.36), for the non-stochastic and 
stochastic parts of the bias, for completely general non-Gaussian initial conditions parametrized 
by the A'-point cumulants ^^^^(fci, ■ ■ ■ ,kj^). 

The outline of the paper is as follows: We will begin, in Section 2, with a qualitative explana- 
tion of scale-dependent stochastic bias. In Section 3, we will show how our intuitive understanding 



^No matter how the fluctuations were created, the parameters have to satisfy the Suyama-Yamaguchi inequahty 
> (|/nl)^ [10] (see also [11-15]). This is easy to understand: we can think of /nl as a measure of the large-scale 
correlation between the potential <& and the locafly measured smafl-scale power, /nl ^ ("l>f$J)/("I"^)($3). On the 
other hand, Tnl is a measure of the large-scale variance in the small-scale power, Tnl ~ {^a^'i) c/ {^e) i^s)'^ ■ The 
inequality r nl > (f /nl)' then arises simply as the condition that the correlation coefficient between the small-scale 
power and <I> must be between —1 and 1. 

'We should note that in this paper we are interested in large-scale stochastic bias. On small scales, non-linear 
evolution and astrophysical processes can create local stochasticity, which is not relevant in our study. 

''More precisely, ki is fixed to the large scale k where we are computing the bias, and k2, ■ ■ ■, kN are integrated 
over a broad range of scales near the halo collapse scale fch ~ M~^^^ . 
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is borne out in the barrier crossing model of structure formation. In Section 4, we will illustrate 
these results with explicit examples. In each case, we also derive our predictions in the peak- 
background split formalism. In a companion paper [1], we prove the mathematical equivalence 
of barrier crossing and peak-background split. We present our conclusions in Section 5. Finally, 
Appendix A discusses the convergence of the Edgeworth expansion for local non-Gaussianity. 

2 Stochastic Bias 

Galaxies reside in dark matter halos. For Gaussian initial conditions and at long wavelengths, 
the fluctuations in the density of halos (^h can be expressed as an expansion in the linear matter 
density field 5. At linear order, the two are simply related by a numerical factor — the bias hg — 
i.e. 5h = bg5. This simple bias relation gets modified for non-Gaussian initial conditions, due to 
a coupling between short and long- wavelength modes. The short modes determine the collapse 
of dark matter halos, while long modes modulate the density on large scales, effectively raising 
or lowering the threshold for the formation of collapsed objects. A non-zero three-point function 
affects the variance of the short modes, leading to a dependence of the number density of halos 
on the amplitude of the long modes. For local non-Gaussianity^ this leads to a dependence of 
the halo density on the long-wavelength gravitational potential $ rather than the matter density 
5 oc V^<I>. This leads to a characteristic scale-dependence in the bias relation, A6 oc k~'^ [7]. It is 
this scale-dependence that allows us to trust the large-scale bias as a probe of initial conditions. 
Crucially, the dependence of the halo density on $ is not something that could be mimicked 
by local dynamics. Dynamical processes don't care about the local value of the potential, but 
are only sensitive to tidal forces which are proportional to V^<I> and 6 (essentially this is a 
consequence of the equivalence principle). Any dependence of the small-scale power on $ itself 
can therefore only come from the initial conditions. This is what makes scale-dependent bias such 
a promising probe of early universe physics, despite all the astrophysical uncertainties associated 
with galaxy formation. 

Stochastic bias arises whenever the density of halos is not 100% correlated with the potential 
$ or its derivatives. In order to develop some intuition, we now give a schematic derivation of the 
effect. In the next section, we will upgrade this to a more formal analysis in the barrier crossing 
approach. If we assume that the primordial perturbations are adiabatic, then the formation 
of halos can only depend on local physics of the fluctuations. Nevertheless, long-wavelength 
variations of the number of halos may depend, not only on the local value of the linear density 
field, but on all of its local correlation functions. Assuming only locality, we may therefore write 
the local halo number density as 



where [...] denotes an average over a small region of characteristic size I that is centered around x. 
Long- wavelength fluctuations in the number of halos can then be understood as a Taylor expan- 
sion. 



^In real space, local non-Gaussianity is parametrized as ^{x) — (j>{x) + ft^i^{(f)'^ (x) — (0^)), where (j> is Gaussian. 
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where bg is the Gaussian bias and 



(9 In Tih 



(2.3) 



It is easy to see (e.g. by sphtting all fields into long and short modes), that for local non- 
Gaussianity the short-scale power is modulated by the gravitational potential, [5^] w (1 + 
4/NL$(ic)). This is the origin of scale-dependent bias in local non-Gaussianity. 

Using the expansion (2.2), we can also evaluate correlation functions between two spatially 
separated points x and x' . We use a prime to indicate that fields are evaluated at x' , while fields 
without a prime are evaluated at x. The matter-halo correlation, in a large region of size L ^ i, 
then is 



{55') 



{55') 



while the halo-halo correlation is 

{SM 



mm 

{55') ^3 ' "'"^'^ {55') {55') 



+ 



(2.4) 



(2.5) 



This leads to the possibility that the bias inferred from ((^h*^) is not equal to the bias inferred from 
('^h'^h')- We characterize this so-called stochasticity of the halo bias by the following parameter^ 



r = 



{W) f{5hS') 



{55') 



Using eqs. (2.4) and (2.5), we find 



13' 



2]/\ 



{55') 



{55') 



{55') 



(2.6) 



+ 



(2.7) 



This simple argument gives reliable intuition for the origin of stochasticity. Specifically, we see 
that if a local variation in the amplitude of [5^] (a;) is uncorrelated with 5{x'), then there is no 
extra contribution to the bias in eq. (2.4). Nevertheless, the halo-halo correlation function in 
eq. (2.5) can still be modified by long-wavelength variations in [(5^](a;). Moreover, the result (2.7) 
makes it clear that stochasticity arises from a non-trivial four-point function of the primordial 
potential. In fact, the real space correlation function ([5^](a;)[(5^](a;')) relates to the collapsed 
limit of the four-point function in Fourier space, i.e. lim|fc^_|_fc2|->>o(*^fci*I*fc2^fe3^fc4)- 



3 Predictions from Barrier Crossing 

In this section, we give a formal derivation of stochastic bias using the classic barrier crossing 
method of Press and Schechter [21]. Our goal is to obtain an expression for the stochasticity 
coefficient (2.6) in terms of the cumulants of the smoothed density field. These in turn can be 
related to A^-point functions of the primordial potential and hence contain information about the 
initial conditions. 

^In practice, we also have to subtract shot noise contributions from (5h5h') and {3^5') — see §3.3. 
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3.1 Definitions and Notation 

We begin with some basic definitions and a description of our notation. Let S{x,z) denote the 
hnear density field (to be distinguished by the hat from the non- linear density field 5). The 
linearized Poisson equation relates 6 to the primordial potential ^, 

5k{z) = aik,z)'^>k , (3.1) 

where 

a{k,z)^l-^T{k)D{z) . (3.2) 

Here, T{k) is the matter transfer function normalized such that T(k) — )• 1 as /c — t- and D{z) is 
the linear growth factor (as function of redshift z), normalized so that D{z) = (1 + z)^^ in matter 
domination. For notational simplicity, we will from now on suppress the redshift argument from 
all quantities. We use Snii^) for the linear field smoothed with a top-hat window function with 
radius*^ Rm = (3M/47rpm)i/3, so that 

6,,{x) = [ e-"'-''W,,{k)5k = [ e-^'=-aM(A:)$fc , (3.3) 

Jk Jk 

Where /,(.)^/(0(.), 

w n \ 1 sin(A;flM) - kRu cos(A;i?M) 

" ' {kR;:;^ ' ^ ^ ^ 

and {k) = WM{k)a{k). Let = {5m)'^'^ be the rms amplitude of the smoothed density field, 
and Kn(M) be its n-th non-Gaussian cumulant, 

K„(M) = Mi , (3.5) 



where the subscript 'c' indicates the use of a connected correlation function. Since 6m and ctm 
are defined via linear theory, Kn{M) is independent of redshift. Similar definitions apply to the 
unsmoothed field 5, in which case we denote the variance and cumulants by a and 

Ultimately, we will be interested in two-point clustering statistics. Let x and x' be two points 
separated by a distance r = |a; — a;'|. Moreover, let a prime indicate that the field is evaluated 
at x' , e.g. = 5m{x'). Fields without a prime are evaluated at x. The joint cumulants are then 
defined by 

K™,„(r,M) ^ ^^"^Sf^'P" , (3.6) 



K^Ar,M,M) ^ ^M^%^. (3.7) 



These cumulants can be related to A^-point functions of the gravitational potential, 

(^fci^fc, •••^fe^)c = (2vr)35D(fci2...Jv)4^^(fci,fc2,...,fcjv) , (3.^ 

where fei2...Af = fci + ^2 + • • • + ^tv. 

^The smoothing scale Rm corresponds to the comoving size of halos of mass M in Lagrangian space. 
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3.2 Edgeworth Expansions 

The probability density functions (PDFs) of weakly non-Gaussian random variables have well- 
defined Edgeworth expansions (for a review see e.g. [22]). Consider first the variables 5m and 5'^- 
It will be convenient to define the rescaled fields 

iy=^ and u' = ^ , (3.9) 

with (z^) = (z^') = and (i/^) = {{v')"^) = 1. The cumulants in eqs. (3.5) and (3.7) then become 
i^n = {i^^)c and Km,n = {i^^ {v')^) c- The Edgeworth expansion for the marginal PDF is 

( / -|^^n Qn \ 1 1 2 

^ P3(z^) , where pg{v) = -—e~ 'i'' . (3.10) 

n>3 ^ ) ^ 

The first few terms can be written as 

V{y) = (l + |if3(^) + |i?4(z^) + • • • ) V,iy) , (3.11) 
where the functions Hn{y) are Hermite polynomials 

1 2 f/" 12 

//„(..) = (-1)^6^^^ — e"'^ . (3.12) 
Similarly, the Edgeworth expansion for the joint PDF is 

-1-15^+ E -^^^-■" a.>ng(,.)n • (3.13) 

m+n>3 ^ ^ / 

In Appendix A, we discuss the convergence properties of this expansion. In the next section, we 
will use it to compute halo-halo correlations. 

The matter-halo case is completely analogous: to construct the joint PDF of the variables 6 
and 5[i, we define rescaled variables O = bjo and v' = b'^joM- The joint PDF p(i>, i^') is then 
given by the Edgeworth series (3.13) with the cumulant Km,n replaced by Km,n. 

3.3 Barrier Crossing 

In the simplest version of the barrier crossing formalism [21], halos of mass > M are identified with 
regions where the linearly evolved smoothed density field exceeds a constant threshold value be 
for collapse. The halo number density n\Jyx) is then given by 

nh(^) = 2^e(5„(a;)-5e) , (3.14) 

with G the Heaviside step function. It has been shown numerically that be ~ 1.42 produces good 
results [23], but for our analytical calculations we don't need to specify a particular value for be- 
The fraction of space occupied by regions above the collapse threshold is 

/>oo 

/(M) = / {M v(y) , (3.15) 
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where fc(M) = 5c/ (J m- Using the Edgeworth expansion (3.11), we find'^ 



1 „ 



/(M) = -erfc(^-^j+p,(i/,) 



^3(M) K,{M) 

H2{Vc) H + 



(3.16) 



When interpreting calculations in the barrier crossing model, it must be kept in mind that 
eq. (3.14) for n'^[x) is the number density of halos in Lagrangian space. Our convention through- 
out this paper is that the power spectra -Pmh(^) and Phh{k) are always computed in Lagrangian 
space. In particular, bg denotes the Lagrangian bias. The relevant quantity to compare to 
observations or simulations is the Eulerian bias which, to lowest order, is given by 6^ = 1 + bg. 

The barrier crossing model also neglects shot noise contributions which arise from the finite 
halo number density nh. Throughout this paper, Phh always denotes the halo-halo power spec- 
trum after subtracting the shot noise contribution l/rih. (There are also shot noise, or one- halo, 
contributions to the matter- halo power spectrum Pmh, which are usually negligible, but are a 
leading source of stochastic bias in the Gaussian case [24, 25].) 

3.3.1 Matter-Halo Correlations 

The correlation between the halo field at x' and the dark matter field at x is given by 

M f'°° - f'°° 

ar,M) = — / [d5] / [d6'J6{x)ni,{x')p{5,6[,) . (3.17) 



2p 

In the rescaled variables u = 6{x)/a and = 6m{'3c') / (Tm-, this becomes 

/OO f'OO 
[du] / [du']up{v,u') . (3.18) 
OO J Vc 

It will be convenient to work in momentum space via ^{k,M) = f d^r e*'°'''^(r, M). To describe 
the correlations of halos in the mass bin [M,M + dM], we take derivatives with respect to M. 
The matter-halo power spectrum is then given by 

^'n.h(fc,M) = ^^ J . (3.19) 

To compute the correlation function (3.18), we substitute the Edgeworth expansion (3.13) for 
p[i',v'). Only terms with exactly one z)-derivative survive the integration, and we therefore find 



i{k,M) = apg{i 



(3.20) 



where * denotes a convolution. We see that the matter-halo correlations, or equivalently the 
non-stochastic part of the halo bias, only depend on the following cumulants 



a,M?j a(fe)($fc$9,---<l>gjc . (3.21) 



''in our notation the halo mass function is dui^/dM = — pui/ M{df / dM). 
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Moreover, we note that the large-scale limit, linik^Q ^{k, M), is determined by the squeezed 
limit of the primordial (n + l)-point function [26], 



iim4"^^^(fc,qi,---,qn) • 



fc^O 



(3.22) 



The explicit form of the cumulants depends on the type of non-Gaussianity. We compute 

some examples in Section 4. 

Keeping only linear terms^ in eq. (3.20), we get 



— + E 



n=2 



fln + -- 



where we defined 



KlJk,M) 



(3.23) 



(3.24) 



It was convenient to factor out the Gaussian term kj ^aaM, since at long wavelengths it becomes 
the matter power spectrum 



Ki^^aaM= / a{k)aM{q){'^k'^q) = WM{k)Pmm{k) Pmm{k) 

Jq 

Evaluating eq. (3.19), we find 



(3.25) 



where 



-fmh(^) ^ -fmm(^) 



1 -1 



'» + E(ft.+^»^)/i,» + 



/3n 



(3.26) 



and (3n 



Hn-i{yc) 



(3.27) 



The ellipses in eq. (3.26) stand for terms that are non-linear in the cumulants. For local non- 
Gaussianity, the derivative terms dfi ^/dlnaM will be negligible, but in principle, we can keep 
them (and sometimes we have to). 

Our expression (3.26) agrees with the general formula for the non-stochastic bias given in [26]; 
however, ref. [26] implicitly found that non-Gaussianity cannot generate large-scale stochastic 
bias. In the next section, we will find the opposite conclusion. The disagreement is easy to 
understand: Ref. [26] claims after their eq. (40) that contributions to Phh from cumulants Km,n 
with m,n > 2 must approach a constant as — ^ 0. This is not true for general non-Gaussian 
initial conditions and exceptions to that statement are precisely what causes the effects discuss 
in this paper. 



In Appendix A, we explain that the lowest order cumulants usually dominate and that products of cumulants 
are suppressed. 
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3.3.2 Halo-Halo Correlations 

Next, we consider the correlation between the halo fields at x and x' , 

MM 



^(r, M, M) 



^rm J — oo 



[d6M] / [d6'j^j]nh{x)nh{x')p{6M,S'j^i 



In the rescaled variables v = 5M{x)/aM and v' = 6M{x')/aM, this becomes 



^(r, M, M) 



J poo 

[du] / [dz.V(^,^^')- 



(3.28) 



(3.29) 



Notice that, in principle, we have allowed for two distinct mass thresholds, M and M. The power 
spectrum of halos in the mass bins [M, M + dM] and [M, M + dM] then is 



^ ^ £ak,M,M) f df{M)df{M) \ 
' dMdM \ dM dM J 



(3.30) 



For simplicity, we will restrict the following presentation to correlations of equal mass halos, 
M = M. The power spectrum for a narrow mass bin around M is then given by 



Phh(fc) 



d^^{k,M, M) 



dMdM 



(df{M)\ 
M=M \ dM ) 



(3.31) 



To compute the correlation function (3.29), we substitute the Edgeworth expansion (3.13) for 



^{k,M,M) =pg{u,)pg{u,) 



+ ^{k2,iHi{j^c) + Ki^2Hi{Uc)) + ^ (K3,ii?2(i^c) + l^l,3H2{l^c)) 



(3.32) 



The form of higher-order cumulants, such as ki^2, «^i,3 and K2,2, again depends on the type of 
non-Gaussianity. We compute some examples in Section 4. 

Keeping only the terms linear in Km,n (this approximation will be justified in Appendix A) in 
eq. (3.32), we find 



^{k,M,M) = Pg{Vc)Pg{Vc){Kl^iaM(yM) 



n=2 

oo oo 



1 Hn 

Jl,n H j Jn,l 



where 



h,n{k,M,M) 

fm,n{k,M,M) 



m=2 n=2 



Kl^n{k,M,M) 
K1,i(/c,M,M)(7m 
Km.a{k,M, M) 

Ki^i{k,M,M)aMcrii 



m 



n 



~ fm,n ~l~ 



(3.33) 



for n > 1 , 
for m, n > 2 



(3.34) 
(3.35) 
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We again factored out the Gaussian contribution, ki^i a^an 



> Pmm{k). Note that /i,^ = fi 



in the large scale limit k <^ ^a/> where /j ^ was defined in eq. (3.24). Substituting (3.33) into 
(3.31), we get 



Phh{k) Pmm{k) 



6J + 26,J:(,3„+^„^ 

d 



n=2 



fl,n 



m=2 n=2 ^ " ^ 



d 



d In cTjy 



fm,n ~l~ 



(3.36) 



Note that, while in the end we always take M = M in this paper, M and M are independent 
variables when calculating partial derivatives of fm,n{M, M). 

3.4 Stochastic Halo Bias 

We now combine the above results to evaluate the stochasticity coefficient 

2 



p 



hh 



Prr 



P 



mh 



Pn 



(3.37) 



where, as usual, it is understood that shot noise is subtracted from Phh and Pmh - Substituting 
eqs. (3.26) and (3.36), we find 



r > 



oo CO / Pi \ / 

m=2n=2 ^ " / \ 



Pn + Pr. 



d 



)fm,n / V ( "I" o 1 ) 



fn,l 



(3.38) 

We note that cumulants Km,n{k) with m,n >2 contribute to the halo-halo power spectrum but 
not the matter-halo power spectrum (3.26), so stochastic halo bias is sourced by these cumulants. 
These cumulants can be written in terms of the (m -|- n)-point functions of the gravitational 
potential. 



Kn,,n{k,M,M) ^ 



M " M 



^ m— 1 „ 

\i=i -^1^ 




(3.39) 



where q = YllLi Qi and q' = Yl^=i 9j- We see that, in general, large-scale stochastic bias 
arises whenever an (m -|- n)-point function C^"^^"^(fci, • • •,fcm+n) is boosted in the collapsed limit 
I Z^i^i ki\ — ?• 0, relative to the product of the corresponding squeezed (m -|- l)-point and (n -|- 1)- 
point functions. In the next section, we will compute eq. (3.38) for a few interesting examples. 



In most cases, we will get stochastic bias from the case m 

,(4) 



n = 2, i.e. a collapsed four-point 
function lim|fc^_|_fc2|->o '^$''(^1) ^2, ^3, ^4) which is larger than the square of the squeezed three- 

(3), 



point function lim^j^o i^i, ^2, ^3) 
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4 Examples 

In this section, we discuss several physical mechanisms that lead to stochastic halo bias. For each 
example, we will derive the result in two different ways: 

1) using a peak-background split (PBS) method; 

2) using the barrier crossing analysis of the previous section. 

We demonstrate explicitly that both approaches lead to the same answers. 
4.1 Tnl Cosmology 

A simple phenomenological way to get a boosted collapsed limit for the four-point function is the 
following generalization of the local ansatz to multiple fields 

$ = Ai(j)i + Bij - {(t)^(|)j)) , (4.1) 

with the Einstein summation convention understood. This structure arises, for example, in the 
curvaton model of [18] (see also [25]), 

$ = + ^ + /^^(l + n)2(V;2-(V;2)) , where ^ =Ii ■ (4.2) 

Here, (p and ip are uncorrelated Gaussian random fields with power spectra that are proportional 
to each other. The three- and four-point functions take the local form 

£!i\kiM, fcs) = /nl [P1P2 + Sperms.] + 0(/i) , (4.3) 
£!i^{kiMMM) = 2(1)' rN4PiP2Pi3 + 23 perms.] +0(r2j , (4.4) 

where we have defined Pi = P^{ki) and Pij = P^{\ki + kj\). However, unlike the single- field 
local ansatz, now Tnl need not be equal to (|/nl) • Instead, the ansatz (4.2) implies t^l = 
(|/nl) (1 + n), in agreement with the Suyama-Yamaguchi inequality, t^l > (|/nl) [10] (see 
also [11-15]). The following limits will be useful in computing the cumulants required in the 



barrier crossing calculation: 

hm (fci, ^2, fca) = 4.UP1P2 , (4.5) 

fcl— !>0 

lim e?^(fci, fc2, fc3, fc4) = 8 (I)' Tnl Pi [P2P3 + P2P4 + P3P4] , (4.6) 

fci— >0 

lim C^^\ki,k2,k;,ki) = 16(|)'tnlPi2PiP3 • (4.7) 



However, before we discuss the explicit barrier crossing result, we present an alternative derivation 
using the peak-background split approach. 
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4.1.1 Peak-Background Split 

PBS is a heuristic procedure for predicting the large-scale clustering statistics of dark matter 
halos. All fields are split into long and short modes — i.e. the Gaussian fields in eq. (4.2) are 
written as (j) = (t)s + 4>t and V' = V^s + V'^- The short scales (< Rm ^ 10 Mpc//i) determine halo 
formation, while the long scales (> 100 Mpc//i) are the ones on which we want to measure the 
clustering of halos. Long modes are therefore always much larger than the Lagrangian size of 
the halos that we consider, i.e. Rm^i <^ 1. The precise split into long and short modes isn't 
important for physical observables, as long as it satisfies the above constraints. 

The long- wavelength modes alter the statistical properties of the small-scale fluctuations. For 
instance, to lowest order, the locally measured small-scale power is o" = cj [1 -|- 2/nl(1 + n)'^^], 
and the locally measured halo number density is 

nk{x) = ni,{6c-5e;a[l + 2f^^{l + U)ije]) ■ (4.8) 
Taylor expanding this expression, we get 

6^=^-^ =bg6e + (3f{l + n)Ui/je , (4.9) 

51nnh , ^ dlnnt , 
5,. and /3/ = 2— . (4.10) 



where 



Hence, we find 



and 



mh 



''j + /3/^)^'mm. (4.11) 



This leads to large-scale halo stochasticity of the form 

-=((i)^N.-/^.)^. (4.13) 



As n — )• 0, this reduces to the classic /nl model, with t^l = (|/nl) and hence no stochasticity. 
4.1.2 Barrier Crossing 

Next, we show that eq. (4.13) can be reproduced precisely from the barrier crossing analysis 
of the previous section. In Appendix A, we show that only the lowest-order cumulants will be 
significant. Here, we calculate the relevant cumulants explicitly: Using eq. (4.5), we get 

= ^4^. (4.14, 

The order-of-magnitude estimates in Appendix A suggest that this will be the dominant contri- 
bution. In particular, we expect, /i^2 > /i,3- We can confirm this explicitly. Using eq. (4.6), we 
get 

fc^O UJ '^NL {/nl = 1) (. . 

/l,3-/i,3 > ^^^-^^ ' 
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where 

J/NL=i)(^^) = \ I I a,,{qi>M{q2)aM{qi2)P^{qi)P^{q2) • (4.16) 



0-3 



Ql 12 

Since 4^^^-'^ is of order A$, we see that the condition /1.3 ^ /i^2 is equivalent to /NL(l+n)A$ <C 
1. This latter condition is always satisfied if all fields are weakly coupled.^ 
Finally, using eq. (4.7), we get 

/2,2^16^i^. (4.18) 
Substituting the above into eq. (3.38) gives 

-=((i)^--/^,)^, (4.19) 

where we have used the relation /3f = 4/^2 = 2(fc — 1), which can be derived by evaluating the 
derivative f3f = 2(91nnh/51n(T in the barrier crossing model [27]. Comparing with eq. (4.13), we 
find that barrier crossing and peak-background split give consistent answers. 

4.2 5rNL Cosmology 

As our next example, we consider a cubic form of local non-Gaussianity.^'^ In this case, the 
non-Gaussian potential is parametrized by the expansion 

^ = + 5nl((/''-3(</.2)0) ^ (420) 

The power spectrum of the non-Gaussian field is 

P^{k) = P^{k) + gl^P^3(k) , (4.21) 

where 

P^^k) =6 [ [ P^{qi)P^{q2)P^{\k - Ql - 92!) • (4.22) 

We note that for scale-invariant initial conditions, (/c3/27r^)P^(/c) = A^, the power spectrum P^z 
is infrared divergent. If the IR divergence is regulated by putting the fields in a finite box with 
length L, then the power spectrum diverges as 

P^z{k) ~ lS/^*^\n^{kL)P^{k) . (4.23) 



^In more detail, to show that /nl(1 + 11) A$ ^ 1, we argue as follows. Assuming that the field ip is not strongly 
coupled, the dimensionless non-Gaussianity parameter /^l' — /nl(1 + Tl)^^^A^ must be < 1. Therefore 

/NL(i + n)A4, = [/nlA*] • [/NL(i + n)3/'A*]'^'<[io-Y'''-[i] = io"' , (4.i7) 

where the bound on the first factor is the current observational bound /nl < 10^ . 

^''We should say from the outset that the large-scale stochasticity in the g^L model will be too small to be 
observationally relevant. Although the non-stochastic and stochastic contributions to Phh(fc) will turn out to be 
parametrically identical (~ g^i^A%P,s>{k)), the non-stochastic contribution is typically larger by a constant factor 
~ 10"*. Nevertheless, the (Jnl example provides an interesting check of our formalism. 
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On large scales, the matter power spectrum therefore is 

Pnnnik) ~ a\k)P^{k) = Pg{k) (l + ISff^^Aj hv^{kL)) , (4.24) 

where we defined Pg{k) = a'^ (k) P(j,{k) . Current observational constraints imply that I^TnlA^I <^ 1. 
To obtain answers to zeroth or first order in ^^nlA^, it suffices to set Pmm — Pg- 

For the barrier crossing analysis, we require the following higher-order correlation functions 

4^)[*'"^^l(fci,fc2,fc3,fc4) = 5NL[AP2P3 + 23perms.] , (4.25) 

4^)ti°°Pl(fc^^ ks, fc4) = [PiP2P^2{ki3) + 11 perms.] , (4.26) 
e^)(fci,fc2,fc3,fe4,fc5,fc6) = 365^ [^1^2^3^^125 + 89 perms.] , (4.27) 

where % = \ki + kj\, Pi = P^{ki), Pijk = P,/,(|fcj + kj + kk\), and 

P^2{k) ^ 2 [ P4q)P4\k - q\) ~ AAlln{kL)P4k) . (4.28) 

Jq 

Note that odd-point correlation functions are zero due to the $ — >• — $ symmetry. Next, 

we will derive the stochastic halo bias both in peak-background split and in barrier crossing. 

4.2.1 Peak-Background Split 

The PBS analysis proceeds as before. Splitting the Gaussian potential into long and short modes, 
4> = (j)£ + (j)s, we find that the locally measured small-scale power is a = cj [l + S^Inl (0£ — i^I))] • 
Moreover, the locally measured value of /nl is /^l = 3(7nl0£ [28]. The halo number density 
therefore is 

n^{x) = nh (<5, -6e;a[l + S^nl - (</>?))] ; ff^) , (4.29) 
where 6i ~ a{kg)(j)£. Taylor expanding this expression, we find 

6h = bg6i + {(pj - {(pj)) + l3ggNL<Pi , (4.30) 

where bg and /3j are the same as in (4.10), and 

ft ^ 3^ . (4.31) 



It follows that 



and 



-Pmh = bgPmm + (3g9Nl.Pm(l, = (^g + 1^9^^^ ' (4.32) 

Phh =[^9+1^9^) Pmm + ^ Pjg^P^^ . (4.33) 



This implies a large-scale halo stochasticity of the form 



-jjaL^ ■ (4.34) 

^ -'mm 
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4.2.2 Barrier Crossing 

We now show that the same result is obtained from barrier crossing. In Appendix A, we argue 
that only the first few cumulants need to be taken into account. It is straightforward to compute 
them explicitly. From eqs. (4.25) and (4.27), we get 

/i,3(A;) ^ ^4^^^='^ and hdk) = [fiMk)? . (4-35) 

This only contributes to the non-stochastic bias. However, since /i 2 = 0) stochastic bias arises 
from /2,2- First, we note that the tree-level four-point function (4.25) leads to a very small and 
scale-independent contribution to /2,2 : 

fiT\k,M,M) ^ (±- [ alMPUl) + A- I o?,M)Pl{<l)] . (4.36) 

Plugging into eq. (3.38) and noting that /32 = and /32 = |, we get a small scale-dependent 
contribution to the large-scale stochastic bias 

(/3| + 2/3/) [ al{q)P'4q)) • (4.37) 



In practice, this contribution to the large-scale stochasticity can't be used as a probe of initial 
conditions, since a contribution to r with r (x \/k (or equivalently a contribution to i-hh(^) 
which approaches a constant as /c — )• 0) is degenerate with other sources of stochasticity such as 
second-order Gaussian bias. Finally, the one-loop four-point function (4.26) leads to the following 
contribution to /2,2 : 

/Pr'(fc) ^ 36,L-#^. (4.38) 

-'mini,"'/ 

The corresponding stochasticity parameter is 

9 2 2 PA^) 1 tA -,0^ 

r[Xoo^\ = lf^f9^^jr^ ^ ' (4-39) 

in agreement with the PBS predictions (4.34). 
4.3 Quasi-Single-Field Inflation 

Our last example is quasi-single field inflation (QSFI) [19]. These models involve extra massive 
scalar degrees of freedom during inflation. In the simplest examples, a single scalar field a of 
mass^^ < mixes with the fluctuation of the inflaton^^ Scf). The mixing communicates 
non-Gaussianity from the hidden (isocurvature) sector to the observable (adiabatic) sector. As 
we now show, it also leads to a significant stochasticity in the halo bias. 



We note that extra scalars with masses close to the Hubble scale H are a natural prediction of supersymmetric 
theories of inflation (see [20] for further discussion). 

^^Recall that S(j> in spatially flat gauge is proportional to the curvature perturbation, = ^5(j). 
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4.3.1 Boosted Four-Point Function 

Again, we need the squeezed and collapsed limits of the primordial correlation functions^^ [14, 
19, 20]: 



k 

where we defined the parameter 



/J \ A 

lim fc2, fca) = P1P2 , (4.40) 

lim ei'^(fci,fc2,fc3,fc4)= 16(1)' TnlIt:^) ^1^3^12, (4.41) 

i2-5>0 y Klii3 J 



3 /9 m2 ^ 

The non-trivial momentum scaling of eqs. (4.40) and (4.41) is a remarkable signature of extra 
Hubble mass scalars during inflation [19, 20, 29, 30]. Moreover, if the mixing between a and 
(or Q is parametrized by a small dimensionless number e < 1, then 

T^^^e-WUf > (I/nl)' . (4.43) 

The enhancement of t^l arises because the trispectrum is generated by the exchange of the cr-field 
which is only weakly coupled to C- The size of the four-point function ((^^) can be estimated from 
the square of the three-point function (C^f) at horizon crossing. 

The boost of Tnl is the result of the small correlation between the curvature fluctuation and the 
massive field, e <C 1. The precise dependence of /nl and Tnl on the fundamental parameters of 
the QSFI Lagrangian can be found in [14]. 

4.3.2 Barrier Crossing 

In QSFI, the higher-order A^-point functions are suppressed by factors of the power spectrum, 
just as in our previous examples. The dominant contributions to the large-scale structure signal 
therefore arise from the squeezed limit of the three-point function and the collapsed limit of the 
four-point function. The relevant cumulants are 

a a[k) af 

and 



Ki,2(/fc) ^ ^{kR,,)'' , (4.45) 



.^Mk) ^ 16 iD'r^.ik'RM'^ Si(A) S|(A) _ 

a («) erf 



A! 



Here, we have defined 

d^A; 



^ Wl{k,)iksRMr^Pmm{ks) , (4.47) 
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See [20] for an intuitive explanation of the scalings in eqs. (4.40) and (4.41). 
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where the integration variable, kg, is one of the short momenta and Rm is the smoothing scale 
defined by eq. (3.4). By definition, Sm(0) = (7^. In the limit A — )• 0, we recover the results of 
the Tnl model. Therefore, we find 



/i,2 = 4/n 



a{k) al 



and 



2,2 



(s/nl)' 



(4.48) 



din (Jm 




Figure 2: Numerical evaluation of eqs. (4.47) and (4.48). For A > 1.0, the cumulant /j j depends signiGcantly 
on the halo mass scaJe M. This is in contrast to local non-Gaussianity, which corresponds to the limit A — > 0. 

To obtain the large-scale stochasticity, we substitute the cumulants into eq. (3.38), 



/32 + /32 



d 



d In (Tm 



/32+/32 



d 



d In cja 



f' 



2,2 



/32+/32 



d 



dlnaA 



fi,2 



(4.49) 



Because the cumulants depend explicitly on Rm, we have to be concerned that the derivatives 
with respect to ctm may this time not be negligible. Indeed, numerical evaluation of the integral 
shows significant cJm dependence of /| 2 (see fig. 2). Keeping the derivative terms, we get 



5)\ 



k^^ 

a^{k) 



din (Tm 



d \ <Si(A) 



DC 



(i) 



f2 

./NL 



^4-2A 



(4.50) 



The characteristic momentum scaling of eq. (4.50) and the natural boost of t^l makes halo 
stochasticity an interesting probe of quasi-single-field inflation. 
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5 Conclusions 



What was the number of hght degrees of freedom during inflation? And, what were their in- 
teractions? The great virtue of primordial non-Gaussianity is that it is sensitive to these basic 
questions about the physics of inflation. In particular, it is well-known that the squeezed limit of 
the primordial three-point function. 



can only be large if more than one light field was dynamically relevant during inflation [31, 32]. 
Remarkably, this statement is independent of the details of the Lagrangian for the inflaton field 
and its initial conditions. Measurements of the squeezed limit therefore have the potential to 
rule out all models of single-field inflation [31, 32]. Moreover, the precise scaling in the squeezed 
limit is sensitive to the details of the mass spectrum [14, 19], allowing a test of extra Hubble 
mass fields, such as those generically expected in supersymmetric theories [20]. Having a large 
three-point function in the squeezed limit modulates the two-point function of halos and therefore 
leads to scale-dependent bias [7]. In the future, this effect may well be our most sensitive probe 
of the squeezed limit. 

In this paper, we have discussed a stochastic form of scale-dependent halo bias. This effect 
arises if the collapsed limit of the primordial four-point function. 



is larger than the square of the squeezed limit of the three-point function. More generally, 
stochastic bias arises whenever a suitable collapsed limit of an (M + A^)-point function is larger 
than the product of the associated squeezed (M + l)-point and (A^ + l)-point functions, where 
M,N > 2. The key tool for obtaining this result, and a main result of this paper, is a pair of 
formulas, eqs. (3.26) and (3.36), for the matter-halo and halo-halo power spectra in a general 
non-Gaussian model parametrized by the A^-point functions of the primordial potential. 

In non-Gaussian models which generate significant stochastic halo bias, the results of this paper 
are important even at a qualitative level. As a concrete example, it should be possible to measure 
/nl and Tnl independently using stochastic bias. This can be done either by measuring multiple 
tracer populations and directly estimating large-scale stochasticity (which has the advantage of 
eliminating sample variance), or from a single tracer population by measuring Phh(^) and using 
the functional form 



to fit for bg, /nl and t^l independently. Recently, ref. [36] showed that if only non-stochastic bias 
is considered, the leading contribution from Tnl is small (in our language, this corresponds to the 
(^(t'nl) contribution to ^1,3) and it is difficult to separate /nl and Tnl, so stochastic bias has an 
important qualitative effect. As another example, in quasi-single field inflation, the stochastic 
bias is larger than the non-stochastic bias by a large factor (parametrically e~^), leading to 
a similarly large enhancement in signal-to-noise when stochastic bias is considered. We defer 
quantitative forecasts incorporating stochastic bias for future work. 






(5.2) 




(5.3) 
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In general, there is no stochastic bias if only a single field (which may or may not be the 
inflaton) generates the primordial curvature perturbation and its non-Gaussianity [16]. Mea- 
suring stochastic halo bias would therefore teach us about the effective number of degrees of 
freedom that generated the primordial fluctuations and its higher-order correlations. In particu- 
lar, stochasticity is sensitive to what we may call "hidden sector non-Gaussianity", i.e. situations 
in which two fields generate the curvature perturbation, but only one (hidden) field is responsi- 
ble for its non-Gaussianity. In this paper, we have derived this effect for general non-Gaussian 
initial conditions. We have also applied our formalism to a number of explicit examples, such 
as curvaton models [18] and quasi-single field inflation [19]. We have shown that halo bias, in 
principle, gives us information about the soft limits of both the primordial three-point function 
and the four-point function. It is therefore a valuable tool in the quest to uncover the physics 
that created the initial perturbations. 
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A Convergence of the Edgeworth Expansion 



In this appendix, we discuss the convergence properties of the Edgeworth expansion for local 
non-Gaussianity. In particular, we will estimate the relative size of the cumulants Hn,m{k, M) for 
general n and m > in the large-scale limit A: — )• 0. For further discussion see e.g. [33-35]. The 
results in this appendix are used in the main text in several places: to justify the approximation 
that non-linear terms in the Edgeworth expansion are negligible in eqs. (3.23) and (3.33), and to 
justify keeping only certain cumulants in the Tnl model (§4.1.2) and the (7nl model (§4.2.2). 

A.l Tnl Cosmology 

We first consider the Tnl model of §4.1. 

Linear terms. — The leading contribution in the A: — )• limit arises from the following contribution 
to the connected correlation function 



I 1 I 1 I 1 I 1 I 1 I 1 I 1 

^ An,m / (V'l(V'V')2 • • • (V'V')n | (V'V')n+l • • • ('i/'V')n+m-lV'n+m)c dK , (A.l) 



where dK = W - OiMi^i) and (V'V')* denotes an auto-convolution evaluated at fcj. The prime 
on the correlation function denotes that we have dropped an overall momentum conserving delta- 
function. The amplitude of the cumulant is given by 

An,m = Cn,^(l + l[fi^'-+rn-2) ^ ^^^^^ ^ ^j^j a^+rn 2 ^^^^^ 

We arrived at eq. (A.l) by using the definition of <I> in eq. (4.2) and expanding out terms to 
produce a connected correlation function. The numerical factor Cn,m in the amplitude (A. 2) 
arises from the sum over equivalent contractions of the fields. The vertical line in (A.l) separates 
the first n terms from the last m. Each contraction gives a factor of P^, and the contraction 
crossing the vertical line carries momentum fe, giving a factor of P^{k) that can be taken out of 
the integral. The power spectrum P^{k) diverges as A: — )• and gives the largest^'^ contribution 
to 

i^n,m- The remaining integral over dK will typically be dominated by the non-linear scale fcnh 
where k'^iPmm{kni) ~ a^,^(A:ni)A| ~ 1. Therefore, we may estimate the integral using ~ Ag^, 
to get 

^i,m ^ c,,m (1 + n)'"-2 A™-2 . for m > 1 , (A.3) 

a[k) 

^^n,m ^ Cn,n. (1 + n)"+— ^ ^ A^+— ^ • ^^^^ for n and m > 1 . (A.4) 

The factor n\m\ appearing in Cn^m is canceled explicitly in the Edgeworth expansion (3.13), and 
as shown in §4.1.2, the condition /nl(1 + 11) A$ <C 1 is always satisfied. This implies that higher- 
order cumulants are subdominant relative to lower-order ones, and hence the only terms we have 
to keep in the Tnl model are ki i , ki 2 = ^2 i and K2 2 ■ 



^■'Subleading contributions arise when both linear ip terms appear on the same side. In such cases, two contrac- 
tions cross the vertical Hne, and the resulting cumulant is finite in the — > limit. 
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Non-linear terms. — When expanding the exponential in the Edgeworth expansion (3.13) we also 
encounter non-linear terms such as k^^{x). First, we will show that, for n and/or m > 1, these 
terms are suppressed by the near-Gaussianity of the primordial perturbations. We distinguish 
two cases: 

• When n > 1 and m > 1, we take powers of the contributions in (A. 4), to find 
where L is an infrared cutoff. This can be written as 

<™ ~ K„,™-<^[/HL(i + n)A<,](^-')("+'"-')(i + n)-^in^-i(fci) • (A.6) 

Using /nl(1 + n)A$ < 1 and (1 + H) > 1, we see that i^nm '^^ suppressed relative to Kn,m 
for n, m > 1. 

• When n = 1 and m > 1, the situation is slightly different. If we take higher powers of the 
results in (A. 3), we find for P > 1, 

- Ki,™-<-^[/..(i + n)A^](^-^)(--^)(i + n)-^a(fc).%Hi(|l) , (A.7) 

Again, as we increase the power P, the contribution is suppressed. However, there is a clear 
difference between P = 1 and P > 1. Nevertheless, in the limit — t- 0, [Pmraik) / a{k)]^^ oc k 
so that these contributions vanish relative to Ki^^n- 

Next, we consider products of the Gaussian piece, k^^. We find for P > 1, 

-fmm(^cq) aP-1/, \ ( \ o\ 

Here, ki^i receives its largest contribution from the peak of the linear matter power spectrum 
^m(^) — k^Pm{k) which occurs at A: = keq, the scale set by matter-radiation equality. Because 
^m(^eq) < 1 at that scale, the modes are still linear and higher powers of ki^i will be suppressed. 
However, in the limit /c — )• 0, ki^i vanishes, while is finite for P > 1. This gives a small 
constant contribution to the halo power spectrum P^h which is a free parameter in practice (we 
discussed this in the context of the ^(nl model in §4.2.2). 

Finally, we look at terms of the form K^^'^n' m'- bound these contributions by using 

the above estimates with the convolution Kn,m*'^n' m' ~ Iq '^rn,ni\Q\) '^n' m'd^"^!)- ^'^^ n,m > 1, 
the convolution will be dominated by the IR, and we find 



'^n,m*^^',m' ~ '^n,m{k)n'^, .^,{k) ^ ^ ^ . (A. 9) 

mmV"'/ 



For m = m' = 1, the convolution is dominated by physics at the non-linear scale, so we may 
simply multiply (A. 7) and/or (A. 8) to find 

<1*4,1 ~ <i(^)4,l(^)^nl • (A.IO) 

As a result, convolutions of different cumulants will be suppressed by /nl(1 + n)A$ ^ 1. 
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A. 2 5NL Cosmology 

Similar arguments apply to the (7nl model of §4.2. 

Linear terms. — First, we note that Kn,m = 0, unless n + m is even. Moreover, only for both n 
and m odd do we get a scale-dependent tree- level contribution to the cumulant. (In the main 
text, we discuss the important special case K2,2-) Schematically, we can write Kn,m ~ 5nl ^X- At 
tree level, we then find 



m — 1 



-2 -fmm(fc) 

a{k) 



Ki,m ~ 5nl' A^-^ • for m odd , (A.ll) 



~ g^f ' A^+"-^ • for n,m odd and > 1 . (A.12) 



Since current observational constraints imply |(/nlA|,| <C 1, the only tree- level terms that we need 
to keep are ki^i, ki 3 = K3 1 and ^3^3. As we discuss in the main text, there is also an interesting 
loop contribution to ^2,2- 

Non-linear terms. — As in the Tnl model, products of cumulants of the form K-nm ^^^^ be suppressed 
due to the near-Gaussianity of the perturbations. The contributions of higher powers of Hn,m is 
nearly identical in both cases: 



When n > 1 and m > 1, we take powers of the contributions in (A.12), to find 



n+m -1 

i/NL ^< 



snlAI, 



^ ' \nP-\kL). (A13) 



Clearly, if 5nlA| ^ 1, then the higher powers of Kn,m are suppressed (if we assume that 
the log is small). 

When n = 1 and m > 1, we take higher powers of the results in (A.ll), to find for P > 1 



5nl A([j A(j, -Pmm(^nl) 



r a2i(^-1)(=¥^) ^mm(fcnl) I ^ ^ A\ 

Again, we find that P > 1 contributions are suppressed by powers of g^i^/^\ <C 1. As in 
the Tnl model, we find that P = 1 has a different scaling with k from P > 1. 



It should be clear that other cumulants behave in the same way as in the t^l model and will be 



suppressed by factors of 5nlA| 
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